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ABSTRACT 

The  three-dimensional  problem  of  wave-trapping  above  a  submerged  round 
sill  was  first  analyzed  by  Longuet-Higgins  (1967)  on  the  basis  of  a  linear 
shallow-water  theory.  The  large  responses  predicted  by  the  theory  were, 
however,  not  well  borne  out  by  the  experiments  of  Barnard,  Pritchard  and 
Provis  (1981)  and  this  has  motivated  a  more  detailed  study  of  the  problem.  A 
full,  linear  theory  for  both  inviscid  and  weakly  viscous  fluid,  without  any 
shallow-water  assumptions,  is  presented  here.  It  reveals  important 
limitations  on  the  use  of  shallow-water  theory  and  the  reasons  for  them.  In 
particular,  while  the  qualitative  features  of  wave-trapping  are  similar  to 
those  of  shallow-water  theory,  the  nearly-resonant  frequencies  differ 
significantly,  and  since  the  resonances  are  narrow,  the  observed  amplitudes  at 
a  given  frequency  differ  greatly.  The  geometry  is  strongly  indicative  of  long 
waves  and  the  dispersion  relation  appears  quite  consistent  with  that,  but  the 
part  of  the  motion  at  wave  numbers  that  are  not  small  has,  despite  the  small 
amplitude,  a  substantial  effect  on  the  response  to  excitation. 
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SIGNIFICANCE  AND  EXPLANATION 


Observations  of  wave  records  at  Macquarie  Island,  in  the  ocean  south  of 
Australia,  have  shown  augmented  response  at  certain  frequencies.  It  has  been 
established  that  these  observations  can  be  explained  as  due  to  the  resonant 
excitation,  of  trapped  modes,  by  the  surrounding  wave  field.  This  phenomenon 
has  prompted  considerable  interest  in  the  possibility  of  wave  trapping. 
Theories,  based  on  the  approximation  of  shallow  water,  lead  to  predictions  of 
very  sharp  resonances  but  these  have  not  been  observed  in  the  laboratory. 

This  paper  develops  a  full  linear  theory,  including  damping,  for  a  special 
geometry,  and  thus  serves  to  define  the  limitations  in  the  results  obtained  by 
using  shallow-water  theory. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


TRAPPING  OF  WATER  WAVES  ABOVE  A  ROUND  SILL 


Yuriko  Yamamuro 


§  1 .  Introduction 

A  theoretical  study  of  the  trapping  that  results  when  a  train  of  small- 
amplitude,  plane  waves  of  a  fixed  frequency  is  incident  on  a  submerged,  steep¬ 
sided,  round  sill  (Figure  1)  was  made  by  Lonquet-Higgins  (henceforth  referred 
to  as  LH)  in  1967.  His  investigation  was  motivated  by  wave  records  taken  at 
Macquarie  Island  showing  the  occurrence  of  regular  oscillations  of  unusually 
large  amplitudes.  In  view  of  these  observations,  LH  considered  a  simplified 
geometry  in  which  the  island  shelf  was  represented  by  a  round  sill,  with  a 
circumference  of  80  km,  submerged  to  a  depth  of  100  m.  He  based  his 
calculation  on  linear,  inviscid,  shallow-water  theory  and  used  separation  of 
variables  in  cylindrical  coordinates  to  determined  the  expressions  for  the 
surface  displacements  for  each  f  the  two  reqions  of  constant  depths.  Recause 
of  the  deoth-independence  of  the  velocity  field  in  shallow-water  theory,  the 
velocity  components  could  not  be  made  continuous  at  the  sill  edge,  and  two 
approximate  matching  conditions  were  used:  the  continuity  of  surface 
elevation  and  the  continuity  of  the  horizontal  component  of  the  mass  flux. 

LH's  (1967)  analysis  showed  the  existence  of  eigenf renuencies  with  very  small 
imaginary  parts.  A  train  of  plane  waves  with  a  frequency  near  such  an 
eigenf requency  could  theoretically  excite  "nearly-trapped"  modes  over  the 
sill,  and  the  response  at  such  modes  was  determined.  The  largest  responses 
were  found  to  occur  at  the  higher  angular  modes  and  at  smaller  ratios  of  the 
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Figure  1 


depths.  These  calculations  have  been  confirmed  and  extended  by  Summerfield 
(1969),  who  applied  shallow-water  theory  to  a  "shelf-island"  model  consisting 
of  a  steep  sided  round  island  rising  from  the  top  of  a  round  sill  of  larger 
radius.  He  showed  the  eigenfrequencies  for  his  system  to  be  closely  related 
to  those  of  LH's  sill  geometry.  However,  laboratory  observations  by  Pite 
(1977)  and  by  Barnard,  Pritchard  and  Provis  (1981)  (henceforth  denoted  by  BPP) 
did  not  reveal  the  large  responses  predicted  by  the  shallow-water  model.  The 
purpose  of  this  paper  is  to  examine  a  complete  linear  theory  in  an  endeavor  to 
explain  these  discrepancies. 

In  §2,  the  equations  governing  the  full  linear  inviscid  theory  are 
presented.  In  §3,  the  velocity  potential  is  calculatd.  This  is  achieved  by 
representing  it  separately  in  the  regions  above  the  sill  and  outside  the  sill; 
the  solutions  in  the  two  zones  are  then  made  to  satisfy  the  necessary 


conditions  of  analytic  continuation  at  the  sill-edge.  This  leads  to  an 


infinite  set  of  linear  equations,  for  which  a  collocation  method  of  solution 
is  described.  In  addition,  an  iterative  method  is  presented  as  a  check  on  the 
extensive  calculations. 

In  §4,  the  theory  is  applied  to  laboratory  conditions  relating  to  the 
experiments  of  BFP.  The  results  show  that  the  modes  that  decay  away  from  the 
sill-edge  and  are  not  included  in  shallow-water  theory,  make  appreciable 
contributions  to  the  wave  amplitudes  above  the  sill.  To  obtain  a  comparison 
on  oceanographic  scales  between  the  full  theory  and  the  shallow-water  theory, 
the  theoretical  predictions  for  the  special  case  of  LH's  sill  geometry  are 
examined  in  §5.  The  differences  are  particularly  striking  near  the 
frequencies  where  the  full  theory  predicts  large  amplitudes,  and  the  reasons 
for  these  differences  are  discussed.  A  theory  allowing  for  weakly  viscous 
effects  is  given  in  §6.  Numerical  computations  were  carried  out  at  a 
frequency  pertaining  to  the  conditions  considered  by  BPP.  At  that  frequency, 
viscous  effects  were  found  to  have  negligible  influence  on  the  amplitudes. 

For  the  entire  range  of  frequencies  investigated  by  BPP,  moreover,  viscous 
effects  are  estimated  to  make  less  than  1%  difference  in  the  amplitudes,  even 
at  the  frequencies  where  peak  amplitudes  are  predicted.  On  oceanographic 
scales,  however,  where  an  eddy  viscosity  may  be  more  appropriate  than  the 
kinematic  viscosity,  damping  may  be  more  pronounced. 


§2.  General  equations 


The  motion  is  assumed  inviscid,  simple-harmonic  in  time,  and  of 
sufficiently  small  amplitude  for  the  use  of  linear  theory.  It  is  referred  to 
cylindrical  coordinates  r,  9,  z  with  the  origin  at  the  undisturbed  water 
surface  (see  Figure  1),  r  measured  outward  in  units  of  the  sill  radius  a, 
and  z  measured  vertically  upward  in  units  of  the  undisturbed  water  depth 
d  above  the  sill;  the  depth  outside  the  sill  is  denoted  by  D  =  d/6.  Let 

x  denote  the  horizontal  coordinate  in  the  direction  9=0.  A  train  of  plane 

* 

waves,  of  frequency  o  and  amplitude  |  p  |  ,  is  incident  on  the  sill  from 

the  positive  x-axis.  The  surface  displacement,  measured  in  units  of  d,  is 

denoted  by  n  and  the  velocity  potential  <j>(r,9,z,t)  is  measured  in  units  of 
2 

d  a  and  satisfies  Laplace's  equation  namely 


3r 


(2.1) 


together  with  the  following  boundary  conditions:  at  the  free  surface,  z  =  0, 


and  ij>  +g/d(J)=0,r>0,0<9<2-rT  ;  (2.2a) 

tt  z 

on  the  horizontal  portion  of  the  seabed,  z  =  -1  and  0  i  r  <  1,  and 
z  =  -1/5,  r  >  1 , 


z 


0 


(2.2b) 
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and  on  the  vertical  wall  of  the  sill  at  r 


1,  -1/S  <  z  <  1 


=  0 


(2.2c) 


At  large  distances  from  the  sill,  the  wave-field  is  assumed  to  consist  of  the 
incident  wave,  whose  surface  elevation  is  the  real  part  of 


n 


i 


.  .  -i(Kx+ot) 
I  hi  T  e 


together  with  waves  that  either  decay  or  radiate  outwards.  Here,  K  is  the 
positive  real  root  of  the  dispersion  relation 


(KD/a)  tanh  (KD/a)  =  Do2/g 


(2.3) 


Different  representations  of  the  solution  will  be  used  in  the  domain  above  the 
sill,  r<  1,  and  in  that  outside  the  sill,  r  >  1.  These  functions  must  be 
analytic  continuations  of  each  other,  for  which  a  necessary  and  sufficient 
condition  is  the  continuity  of  <J>  and  9<)>/3r. 
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§3.  Form  of  solutions 


Let  the  velocity  potentials  for  r  <  1  and  r  >  1  be  denoted  by  <J> 

s 

and  respectively,  and  be  simple  harmonic  in  time  with  radian  frequency 

o.  Make  the  decomposition  ^  ^  +  <}>R  where  $  represents  the  incident 

plane  wave  train  and  A  represents  waves  generated  by  the  sill.  A  is 

R  X 

calculated  by  separation  of  variables  in  Cartesian  coordinates: 


4> 


X 


-i(  aXQx/D+cjt) 
e 


cosh  X0(6z+1) 
XQ  sinh  XQ 


+  * 


(3.1) 


Here,  and  throughout  the  paper,  the  notation  +  *  will  be  used  to  denote  the 
addition  of  the  complex  conjugate  of  all  preceding  terms.  The  free-surface 
boundary  conditions  require  XQ  to  be  the  real  and  positive  root  of  the 
dispersion  relation  (e.g.  see  Davis  &  Hood  1976) 


X  tanh 


X  =  Do  /g 


(3.2) 


For  convenience,  the  expression  for  <*>  is  converted  to  cylindrical 
coordinates  (Bateman  Manuscript  Project  1953,  Vol.  2,  §7.2.4,  p.  7,  Eqn.  27): 


_ .  cosh  X  ( <5z+1 )  00 

j>  =  e  10  — — - :  -  -  —  V  e  i  mJ  (aX  r/D)cos  m0  +  * 

I  2  X„  sinh  X„  '■  m  m  0 

0  0  m=0 


(3.3) 


where  e  =  1  ^  m  ,  !!  For  the  purpose  of  computations  for  the 

m  2  if  m  j*  0 

linearized  problem,  Inl^  normalized  to  unity. 

The  generated  field  in  the  region  outside  the  sill,  r  >  1,  may  be 
represented  as  the  sum  (Havelock  192^) 


-6- 


-icrt 


*R  -  6 


T  cos  m9 {B  H  (aL r/D)cosh  X.(5z+1) 
mO  m  0  0 

m=0 

(3.4) 

00 

+  T  B  K  (aA  r/D ) cos  X  (6z+1)/K  (aX  /D)}  +  *  . 

mn  m  n  n  m  n 

n=1 


Here,  elements  of  {±iX  :  m  =  1,2,...}  are  the  purely  imaginary  roots  of  the 

m 

dispersion  relation  (3.2).  Note  that  the  model  expansion  is  complete  (Davis  & 

2 

Hood  1976)  and  that  X  =  md  for  m  sufficiently  large  or  for  Da  /q 

m 

sufficiently  small.  The  notations  for  the  Bessel  functions  are  those  of 
Abramowitz  and  Stegun  (1972,  §9).  Each  term  in  (3.4)  satisfies  the 
appropriate  conditions  at  infinity  as  well  as  the  free-surface  conditions  and 
boundary  conditions  on  the  seabed.  The,  as  yet  undetermined,  coefficients 
Bm0  and  B^fn?*  0)  are,  respectively,  coefficients  of  each  radiating  mode 
and  each  non-radiating  mode,  which  is  confined  near  the  sill-edge.  From  now 
on,  the  subscript  on  XQ  and  the  superscript  on  the  Hankel  function  will  be 
dropped. 

The  solution  in  the  sill-reqion,  r  <  1,  may  be  represented  as  the  sum 
of  solutions  obtained  by  separation  of  variables,  namely 


QO 

<(>  =  e  1<Tt  T  cos  m0{A  J  (akr/d)cosh  k(z+1) 

vs  mO  m 

m=0 

OO 

+  ^  A  I  (.!«  r/d)cos  k  (z+1  )/I  (ak  /d)}  +  * 

mn  m  n  n  m  n 

n=  1 


(3.6) 


Here,  elements  of  {k  ,  ik 


1,2,...}  are  the  roots  of 


2 

k  tanh  k  =  da“/g 


(3.6) 


of  which  each  of  the  are  positive.  Each  term  satisfies  the  appropriate 

conditions  at  the  edge  of  the  sill,  at  the  free-surface  and  on  the  seabed. 
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The,  as  yet  undetermined,  coefficients  A^q  and  Amn(n/0)  are,  respectively, 
coefficients  of  each  wavelike  mode  and  each  spatially  decaying  mode.  In  what 


follows,  the  subscript  on  k^  will  be  dropped. 

The  continuity  of  3ij)/3r  at  r  =  1,  -1  <  z  <  0,  0  <  9  <  2ir,  yields  one 
of  the  two  conditions  at  the  sill-edge: 


3r 


34>e 

3r 


for  -1  <  z  <  0 


(3.6] 


0  for  -1/S  <  z  <  -1 


By  virtue  of  the  orthogonality  property  of  the  set  {cosh  X(6z+1), 

cos  X  (6z+1)  :  n  =  1,2,...},  this  condition  yields,  for  m  =  0,1,... 
n 


[B  H'  (aX/D) 
mO  m 


F  ] 
n 


aX 


h(  X) 


3<j>  . 
m 

3r  r=1 


cosh  X( 6z+1 )dz 


(3.7) 


Here  d>  denotes  the  coefficient  of 
m 

h(X)  are  defined  in  the  appendix. 

K'  (aX  /D)  aX  „  3<f>  . 

m  n _ n  _  ,0  m 

ran  x  (aX  /D)  D  (1  n  '  - 1  3r  r= 
m  n 


A  -iot 

cos  m0  e  in 

For  n  =  1,2,..., 

.  cos  X  ( Sz+1 )dz. 
1  n 


* 

A  and  F  and 
s  m 


The  second  condition  at  r  =  1 ,  - 1  <  z  <  0 ,  0  <  9  <  2it,  is  the 

continuity  of  i>.  This  yields  a  second  relation  between  B_  and  A_  ,  from 

'  mn  mn 

which  the  Bmn  are  eliminated  by  means  of  equations  (3.7).  The  resulting 
equations  for  the  Amn  are: 


A  {cosh  k(z+1)J  (ak/d)  -  £  A(z)}  +  ^  A  {cos  k  (z+1)  -  £  (z)} 

mu  m  mO  mp  p  nm 

D=1 

*  H  ( aX/D) 

=  {F  -  F  ■— — - }  cosh  X(6z+1)  for  -1  <  z  <  U,  m  =  0,1,... 
m  m  H  (  a  a/D  ) 
m 


(3.8) 


-  R  - 


The  notation  is  defined  in  the  appendix.  When  no  spatially  decaying  modes  are 

included,  equation  (3.8)  is  first  multiplied  by  cosh  k(z+1)  and  then 

integrated  over  z.  With  the  inclusion  of  the  decaying  modes,  a  collocation 

method  is  used  to  solve  the  resulting  infinite  number  of  equations:  that  is, 

an  N  by  N  matrix  equation  is  constructed  by  the  application  of  equation 

(3.8)  at  N  values  of  z  in  the  range  -1  <  z  <  0  and  by  neglecting  the 

higher-order  decaying  modes  {A  :  p  >  N } . 

mp 

Justification  for  the  matrix  truncation,  which  is  performed  in  the 
calculations  of  the  following  sections,  is  as  follows.  For  the  investigation 
on  the  laboratory  scales  described  in  §4,  the  N  values  of  z  were  chosen  to 
be  {z  :  z  =  -i/(N+1),  i  =  1,...,N}  for  N  =  2,6,11,16,21,41.  It  is  shown  in 

§4  that  the  differences  in  the  resulting  amplitudes  for  N  >  6  was  of 
negligible  importance.  Hence,  it  is  apparent  that  the  above  choice  for  the 
values  of  z,  with  N  =  6,  was  adequate  for  computations  involving  similar 
sill  geometries.  The  computations  in  $6  concern  scales,  for  which  the 
parameter  a/d  was  much  greater  than  in  §4,  s"  -'-at  the  local  modes  decayed 
faster  away  from  the  sill-odue  than  r •:>•••*»  -*  ^4.  nee  experimental  results 

by  Pite  (1977)  and,  independently,  !v  PM  in:;  almost  no  presence  of 

exponential  decay  in  the  wave-fid  t  -,ur  d  .•  Ml  -.>  ige,  ;t  •.  s  expected  that 
the  faster  the  local  modes  decay,  t‘-<-  Mss  will  he  their  response.  Therefore, 
N  =  6,  which  was  sufficiently  accurate  for  ^4,  was  also  expected  to  give 
suitable  accuracy  for  §5,  and  has  also  been  used  for  the  computations 
described  in  ^5. 

The  computations  are  made  difficult  by  the  presence,  in  e Tuition  (1.9), 

of  the  infinite  series  in  the  terms  P.  (z)  for  n  =  0,1,...,  m  =  0,1,..., 
a  mp 

which  are  slowly  convergent.  Because  of  this,  a  numerical  check  on  the 
complicated  comnut a t.  ions  was  required,  and  an  iterative  method  is  now 
presented  for  this  purpose. 


-O- 


The  decomposition  of  <j)g  in  the  radial  variable  yields  (as  shown  by 

equation  (3.5))  a  principal  part  Jm(akr/d),  which  is  wavelike,  and  an 

infinite  number  of  modes  Im(aknr/d)  that  decay  away  from  the  sill-edge.  If 

the  value  of  the  principal  radial  eigenfunction  at  the  sill-edge  is  not  too 

small,  then  an  iterative  procedure  may  be  adopted.  From  now  on,  the 

coefficient  of  cos  m0  e  1CTt  in  <b  and  will  be  denoted  by  A  . 

A  preliminary  outer  flow  is  determined  from  the  radiation  condition  and 

the  boundary  condition  3<f>m/3r  =  0  at  r=1,  -1  <  z  <  0.  This  outer  flow 

then  determines  the  boundary  conditon  at  r  =  1,  -1  <  2  <  0,  for  a  sill  flow, 

through  the  continuity  of  This  sill  flow  then  yields  3<j>  /3r  at  r  =  1, 

m 

-1  <  z  <  0,  from  which  a  second  outer  flow  is  calculated,  and  so  forth. 

Equations  (A)  and  (B)  of  the  appendix  arise  from  this  scheme  and  are  to 

be  used  as  follows.  When  the  Amn  in  ^  are  known,  the  coefficient  Bmn  of 

d>  are  calculated  from  the  equations  (A).  Note  that  here,  H'(a\/D), 

^0  m 

I  (ak  /d)  and  K'(a\  /D )  do  not  vanish.  When  the  13  are  known,  the 
m  n  m  p  ~n 

coefficients  in  A  are  calculated  from  the  equations  (B) .  This  iteration 

s 

diverges  near  the  zeros  of  Jm(ak/d),  which  lie  near  the  peak  periods.  This 
scheme  therefore  fails  for  the  parameters  of  most  interest,  but  it  was  felt 
important  to  use  the  method  to  provide  a  check  on  the  computations  made  with 
(3.8). 


-in- 


§4 .  Application  to  laboratory  scales 

This  section  concerns  some  calculations  made  to  correspond  to  the 
laboratory  scales  used  by  BPP.  For  their  experiments,  a  =  50  cm, 
d  =  1.75  cm,  D  =  15.4  cm  and  the  forcing  periods  ranged  between  0.75  and  1.20 
seconds. 

The  full  linear  theory  with  no  spatially  decaying  modes  was  used  to 
compute  the  response  curve  for  each  model  number  and  were  compared  with  the 
linear  shallow-water  theory.  Figures  4.1(a)  and  4.1(b)  display  the  results 
for  the  0th  mode  and  figure  4.2  those  of  the  6th  mode.  The  curves  show  that 
the  peaks  in  the  lower  freauency  range  occurred  near  those  of  linear  shallow- 
water  theory.  At  small  modal  numbers,  however,  the  response  was  significantly 
diminished  and  at  large  modal  numbers,  the  bandwidths  were  small  so  that  the 
LH  theory  was  not  useful  in  calculations  of  the  amplitudes.  In  the  higher 
frequency  range,  there  was  no  evidence  of  the  peaked  response  predicted  by  the 
LH  theory.  In  this  connection,  note  that  there  are  a  large  number  of  zeros  in 
the  radial  eigenfunction  Jm(akr/d)  in  the  sill-region.  Thus,  the  horizontal 
scale  of  the  motion  is  so  much  smaller  than  the  radius  that  the  motion  is  not 
a  long  wave,  and  this  is  the  most  probable  source  of  the  substantial  error  in 
that  theory. 

Computations,  accounting  also  for  the  decaying  modes,  were  carried  out  at 
the  forcing  period  of  1.181  seconds  where  the  measurements  by  BPP  yielded  the 
largest  amplitude  above  the  sill.  The  difference  between  the  empirical  values 
(In  (r.,  0 . ) I  :  r.  =  1/14  +  (i-1)/7,  i  =  1,...,7,  0.  =  ( i- 1 > / 1 ft , 

Fill  l 

i  =  1,...,1ft)  and  the  theoretical  values  ( | | )  of  the  wave  amplitudes  was 

measured  by 
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14.00 


E2  = 


l  1  WW1  -  lVri'ej,l)2ri 

i=1  j=1  J 


7  19 


I  I  wh  1 WV1  ri 

i=1  j=1  J 


1/2 


(4.1  ) 


and 


max  |lnE,(ri»0j)l  ~  l9T(r^/0j)l 

i/j 


Em  max  I n_(r  .  ,  0  . )  I 
.  .  fc-  1  J 

1.3 


(4.2) 


The  coefficients  w^  and  w‘  .  represent  the  weights  determined  by  Simpsons 


rule  for  integration  over  9  (0  <  0  <  it)  and  t(r^  <  r  <  ry^'  together  with 
additional  terms  from  the  integration  over  0  <  r  <  and  r?  <  r  <  1 .  For 
example,  the  square  of  the  numerator  of  E2  is 

{S(r1)r1  +  S(r7)r7  +  4(S(r2>r2  +  S(r4>r4  +  S(rfe)r6> 


-  2<S(r3)r3  +  S(r5)r5)}/21  +  (Sfr^r,  +  S(r7)r7)/14  , 


(4.3) 


where,  for  i  -  1,...,7, 


S(  r  .  )  =  tt/54  {4 


j=2n , n= 1 , . • .9 


(lnE(x.,0.)1  -  lnT(r.,0j)l)‘ 


+2  l  ( lnE(r.,0.)|  -  lnT(r.,e.)|)' 

j=2n+ 1 , n= 1 , . . .8 


+  <  I  Vri' V  '  1  W  V  1  r 


(4.4) 


(  ,nE(ri'019)l  •  'W 


(r;  ,  ^l9)  |  )*•}  . 


A  consistency  check  of  E2  made  by  using  the  w^  .  according  to 

the  midpoint  rule. 
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When  only  the  travelling  modes  A^fm  =  0,...f8)  were  included  in  the 
computations,  the  improvement  in  Ej  over  the  LH  theory  was  8%.  When  the 
first  decaying  mode  Am.j(m  =  0,...,8)  was  included,  together  with  a  large 
number  of  "outer"  decaying  modes,  E £  was  further  improved  by  13%.  With  the 
inclusion  of  from  5  to  40  "sill"  decaying  modes,  E2  decreased  by 
negligible  amounts.  These  results  document  the  importance  of  the  decaying 
modes,  which  had  previously  been  though  negligible  (Pite  1977). 

The  above  computational  results  are  explained  as  follows.  The 

experimental  conditions  may  be  modeled  by  the  full  equation  (3.8)  when  6  and 
2 

Do  /g  are  both  small.  Under  these  conditions,  the  response  equation  becomes 
approximately: 


H  (aX/D) 
m 


A  {J  (ak/d)  -  k  J'  (ak/d)  ■  -■  +  X  +  Y)} 

mO  m  m  \H  (a X/D ) 

X  m 


(4.5) 


*  H  (aX/D) 
j,  _  _m _ 

‘  m  m  H '  ( a  X/D ) 
m 


where  X  = 


y 

3  r2 

it  6  n=  1 


sin  nw(  1  -  6) 


l 


sin^  VTT (  1  -  5)  ^ 


5  2  2  2 

u  p=1  n=1  n(n  6  -p  ) 


X  is  (1/6  )  and  originates  from  the  "outer"  decaying  modes,  of  which  a 
large  number  must  be  included  in  the  computations.  Y  is  0(1)  and 
originates  from  the  "sill"  decaying  modes.  Equation  (4.5)  indicates  that  the 
inclusion  of  the  "outer"  decaving  modes  is  crucial  to  the  satisfaction  of  the 
sill-edge  conditons,  but  that  the  inclusion  of  the  "sill"  decayinq  modes  is 
not.  This  accounts  for  the  differences  in  Ej  iecribed  in  the  last 
paragraph . 
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The  maximum  amplitudes  above  the  sill  and  the  amount  of  mechanical  energy 
over  the  sill  were  computed  using  the  full  linear  theory  with  5  "sill"- 
decaying  modes.  The  results,  for  the  experimental  range  of  forcing 
frequencies,  are  shown  in  figures  4.3  and  4.4  together  with  those  of  the  LH 
theory.  Compared  with  the  peaks  given  by  the  LH  theory  for  this  range  of 
frequencies,  those  of  the  full  linear  theory  were  shifted  significantly  to 
higher  periods,  with  increased  magnitude  in  the  case  of  the  'energy'  graph. 

The  peaks  of  the  full  linear  theory  occur  near  the  complex  zeros  of  the 
coefficient  of  in  equation  (4.5).  The  LH  theory  peaks  near  the  zeros, 

k  J'  (ak/d)H  <aX/D)  2  2 

k,  of  J  (ak/d)  -  - ^  — re— ; -  =  0  where  1;  =  do  /q  =  5X  .  This 

m  X  H '  (  a  X/D ) 

m 

2 

indicates  that  the  shifts  arise  through  the  terms  X  +  Y  which  is  0(1/5  ) 

2 

for  small  5.  The  LH  theory  is  therefore  useful  here  only  when  k/6  <<  1. 

If  ak/d  is  too  small,  however,  there  is  no  resonance  in  the  system.  This 

2  2 

indicates  that  the  linear  shallow-water  theory  is  useful  if  Do  /g  and  k/5 
are  both  small  and  ak/d  is  not  small. 

It  apnears,  therefore,  that  a  shallow-water  theorv  is  not  necessarily  a 
good  approximation  to  the  full  linear  theory  even  thouqh  the  boundary 
conditions  seem  to  indicate  it.  The  reason  lies  in  the  extreme  sensitivity  of 
the  matrix  equation  arisinq  from  (3.8)  (and  hence  the  response)  to  what  miqht 
he  presumed  to  be  snail  perturbations.  for  example,  the  terms  arisinq  from 
the  decaying  modes  might  at  first  guess  be  thouqh  to  he  unimportant,  indeed, 
ti.eir  relative  excitations  for  P  =  1,2,...  m  =  0,1,...  can  he 

shown  to  be  small.  Also,  the  wave  numbers  k  and  X,  when  evaluated  hv  the 
full  linear  dispersion  relations,  differed  only  sliuhtlv  from  their  shallow- 
water  values.  T|ipir  small  differences,  however,  caused  noticeable  effects  on 


(3.8)  through  the  Bessel  functions,  especially  at  forcing  frequencies  at  which 
the  resulting  matrix  equation  was  nearly  singular,  and  these  are  just  the 
frequencies  of  peak  response.  The  standard,  diagnostic  checks  on  dispersion 
relation  and  modal  amplitudes  do  turn  out  to  mislead. 

An  implication  of  this  result  is  that  shallow-water  theory  should  be  used 
only  with  great  care  in  interpreting  experimental  results  concerning  motions 
that  in  theory  exhibit  sharp  changes  in  the  response  for  small  changes  in  the 
forcing.  Another,  is  that  great  care  may  need  to  be  taken  with  the  numerical 
work  at  such  frequencies. 
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§ 5 .  Application  to  oceanographic  scales 


LH  (1967)  has  calculated  the  response  | An | ,  where 


oo  ^  ^ 

n  =  V  0t*A  J  (k,r  ),  k,  =  a/(gd)^  for  a  selection  of  frequencies. 

u  n  n  i  i 

—  00 


The  corresponding  response  using  the  full  linear  theory  is  |C  |/e  ,  where 

mu  m 


n  =  e 


.  <*>  <•>  I  ( ak  r/d ) 

1  )  cos  mere  J  (akr/d)  +  V  C  n  ■  ]  +  *,  C  .  . 

mO  m  mn  I  (ak  /d)  mO  - 

m=0  n=1  m  n 


2  2 
A  (ido  /g)  cosh  k,  m  =  0,1,...  and  C„  =  A  ( ido  /g)cos  k  ,  m  =  0,1,..., 
mo  mn  mn  n 


n  =  1,2,...  .  The  results  for  the  frequencies  and  depth-ratios,  5, 
considered  by  LH  are  listed  in  Table  5.  For  <$  =  9/16,  there  were  not  very 
large  differences  between  the  predictions  of  the  two  theories.  At  6  =  1/16, 
however,  the  LH  theory  predicted  very  large  responses,  especially  at  larger 
values  of  m.  The  matrix  equation  was  then  sensitive  to  approximations 
inherent  in  shallow-water  theory,  mentioned  toward  the  end  of  §4,  and  the  two 
theories  produced  very  different  results. 

It  is  of  interest  to  delineate  the  parameter  ranges  for  which  the 


shallow-water  theory  may  be  useful.  This  theory  is  based  on  the  smallness  of 
2 


the  parameter  Do  /g  and  the  practical  range  for  which  the  theory  appears  to 
be  useful  is  limited  to  values  of  this  narameter  rather  less  than  0.1  (Table 


4.1,  Silvester  1974).  On  the  other  hand,  the  dimensionless  frequencies 


fa 


ao/(gd)  i  at  which  resonance  occurs  lie  close  to  the  zeros  of  Jm(z);  the 


smallest  zeros  of  increase  as  m  increases  and  for  m  greater  than 


6,  the  first  zeros  are  greater  than  10.  Since  the  shallow-water  parameter 


.V, 


2  I/22 

is  related  to  an/(gd)  ^  by  Dr/Yg  =  (ao/(  gd)  ^  )  (d/a)  /iS,  the  LH  theorv  is 

likely  to  be  useful  for  an  appropriate  combination  of  (i)  d/a  small,  (ii) 

<5  not  too  small,  and  (iii)  m  not  too  large.  The  large,  narrow,  'resonant' 


-2 


of  full  linear  theory  «itb  Lcuuuut-Hiegu.u  rarronK 


| A  I  (1 967)  in  the  surface  elevation 

'  n’ 


Mode  no. 


Dimension! ess  frequency  Longuot-lUgg ms 

art  1  1 


Full  linear  theory’s 
coeff.  of  travelling 
mode 


4 . 921 
8.219 

7.447 

10.798 

14.027 

9.821 

13.408 

16.723 

19.922 

23.235 

12.121 

15.890 

19.352 

22.662 

25.869 

29.083 

32.480 


16.23 

8.057 

91.37 
18.90 
8.  593 

581.1 

85.81 

24.21 

10.42 

7.122 

3895 

437.8 

98.12 

31.76 

13.55 

7.954 

6.512 


11.88 

7.267 

2.145 

5.756 

6.233 

0.2798 

1.176 

2.840 

4.331 

4.196 

0.0757 

0.2235 

0.7472 

1.714 

2.797 

3.269 

2.769 


4.796 

5.632 

6.799 

8.021 


9.240 

12.832 

11.630 

15.025 


2.871 

3.269 

3.546 

4.447 

5.841 

3.738 

10.83 

3.018 


2.775 

2.921 

3.526 

4.482 

5.644 

2.516 

5.351 

3.026 


6.8903 


0.978 


1.357 


2 


peaks  predicted  by  the  LH  theory,  however,  occurred  for  parameter  ranges  not 


satisfying  these  restrictions.  For  instance.  Table  5  (d/a  =  0.008)  shows 
that  the  most  spectacular  peaks  occurred  for  the  larger  m  and  smaller  6 
and  that  these  peaks  were  not  predicted  by  the  more  accurate  linear  theory. 


* 

1 
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§6.  Weakly  viscous  effects 


Consider  a  periodic  flow  above  the  sill  at  a  frequency  a*  At  the  solid 

Vo 

boundary,  there  will  be  a  Stokes  boundary  layer  of  thickness  ( v/a)  ,  where 

v  is  the  kinematic  viscosity  of  the  fluid.  For  laboratory  scales,  v  is 

?  - 1 

assumed  to  be  0.01  cm  sec  (for  water)  and  for  the  field  scales,  v  may 

2  - 1 

be  the  eddy  viscosity,  perhaps  of  the  order  of  10  cm  sec  .  In  the  sill 

Vo 

geometry  of  BPP,  (v/a)  yd  =  0.05  and  for  the  case  considered  by  LH,  this 
parameter  was  an  order  of  magnitude  smaller.  The  effect  of  viscous 
dissipation  at  the  bottom  boundary  outside  the  sill-region  is  expected  to  be 
small  compared  with  that  inside.  In  order  to  estimate  the  effect  of  viscosity 
on  the  total  flow,  it  is  therefore  plausible  to  take  the  boundary  layer  into 
account  for  r  <  1,  and  to  neglect  it  for  r  >  1.  Following  the  work  of 
Mahony  and  Pritchard  (1980)  the  boundary  condition  to  be  posed  at  z  =  -1, 
r  <  1,  0  <  0  <  2tt,  is  found  to  be 


3*  iir/4  3*  .  3  4 

1*  C1  *  7l*0,V2) 

oZ 


v>  2 

where  e  =  ( v/a)  yd  and  =  (d/a)  . 

Computations  for  the  sill  geometry  studied  by  BPP  showed  approximately  a 
1%  difference  in  E 2  (as  defined  by  4.1)  from  the  inviscid  case,  over  the 
entire  frequency  range.  Thus  it  would  appear  that  the  effect  of  viscosity  is 
not  important  under  their  conditions.  The  suggestion  by  Dite  (1977),  who  used 
a  guasi-empirical  theory  which  involves  a  fluctuating  body-force,  that  viscous 
ffects  were  significant  on  his  laboratory  scales,  does  not  appear  to  be 
justifiable  since  his  values  of  and  were  similar  to  those  of  BPP. 
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APPENDIX 
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i  (z) 
mp 


m^d 
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for  p  ^  0 


Equations  (A) 


For  m  =  0,1,..., 


8mf>  =  7777 1  TmT  (tt,  [kA  ng(  X,k)  J'  (ak/d) 
m0  H ' ( aX/D )  3  Xh { X  y  m0  m 


<*>  I*  (ak  /A)  t 

+  Y  k  A  Tm,  g(  X,  ik  )]  -  F  } 

L „  n  inn  I  (ak  /A)  n  n 

n=  1  m  n 


and  for  p  =  1,2,..., 


B 


K  (aX  /D) 

-  ~ r; — -  -  .  .  {kA  g(  i  X  ,k ) J 1 (  ak/d) 

mp  K'(aX  /D)  6X  h( lXp)  m0  pm 

m  p  p 


co  I '  (  ak  /A) 

T  k  A  g(iX  , ik  )} 

L.  n  mn  I  (ak  /A)  ’  " 
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Equations  (B) 


For  m  =  0,1,.,. 


a-  =  tt— rTTrmrr  (ts  .h  (ax/D)  +  f  ]g(X,k) 
j  {ak/d)f(k)  mO  m  m 


I  B  g(  i X  ,k)  } 
L  mn  n 


Here  f(k),  Fm  and  g(X,k)  are  defined  in  the  appendix. 


For  p  =  1,2,..., 


{[B  H  (aX/D)  +  F  ]g(X,ik  ) 
mp  f ( lk  )  mO  m  m  p 


+  y  B  g(iX  ,  ik  )} 
“  mn  n  p 
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The  three-dimensional  problem  of  wave-trapping  above  a  submerged  round  sill 
was  first  analyzed  by  Longuet-Higgins  (1967)  on  the  basis  of  a  linear  shallow- 
water  theory.  The  large  responses  predicted  by  the  theory  were,  however,  not 
well  borne  out  by  the  experiments  of  Barnard,  Pritchard  and  Provis  (1981)  and 
this  has  motivated  a  more  detailed  study  of  the  problem.  A  full,  linear  theory 
for  both  inviscid  and  weakly  viscous  fluid,  without  any  shallow-water  assumptions, 
is  presented  here.  It  reveals  important  limitations  on  the  use  of  shallow-water 
theory  and  the  reasons  for  them.  In  particular,  while  the  qualitative  features  — . 
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of  wave-trapping  are  similar  to  those  of  shallow-water  theory,  the  nearly- 
resonant  frequencies  differ  significantly,  and  since  the  resonances  are  narrow, 
the  observed  amplitudes  at  a  given  frequency  differ  greatly.  The  geometry  is 
strongly  indicative  of  long  waves  and  the  dispersion  relation  appears  quite  con¬ 
sistent  with  that,  but  the  part  of  the  motion  at  wave  numbers  that  are  not  small 
has,  despite  the  small  amplitude,  a  substantial  effect  on  the  response  to 
excitation . 
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